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This study presents a new class of turbulence model designed for wall bounded, high 
Reynolds number flows with separation. The model addresses deficiencies seen in the 
modeling of nonequilibrium turbulent flows. These flows generally have variable adverse 
pressure gradients which cause the turbulent quantities to react at a finite rate to changes 
in the mean flow quantities. This ’’lag” in the response of the turbulent quantities can’t 
be modeled by most standard turbulence models, which are designed to model equilibrium 
turbulent boundary layers. The model presented uses a standard 2-equation model as the 
baseline for turbulent equilibrium calculations, but adds transport equations to account 
directly for non-equilibrium effects in the Reynolds Stress Tensor (RST) that are seen in 
large pressure gradients involving shock waves and separation. Comparisons are made to 
several standard turbulence modeling validation cases, including an incompressible bound- 
ary layer (both neutral and adverse pressure gradients), an incompressible mixing layer 
and a transonic bump flow. In addition, a hypersonic Shock Wave Turbulent Boundary 
Layer Interaction with separation is assessed along with a transonic capsule flow. Results 
show a substantial improvement over the baseline models for transonic separated flows. 
The results are mixed for the SWTBLI flows assessed. Separation predictions are not as 
good as the baseline models, but the over prediction of the peak heat flux downstream of 
the reattachment shock that plagues many models is reduced. 

Greek Symbols 

a Angle of attack or turbulence modeling constant 
/ 3 Clauser’s pressure gradient parameter, turbulence modeling constant 

5 Boundary layer thickness, m 

6* Displacement thickness, m 
e Dissipation per unit mass, ?n 2 /s 3 

k Turbulence modeling constant 

H Dynamic viscosity, kg/m ■ s 

v Kinematic viscosity, m 2 / s 

uj Specific dissipation, s^ 1 
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p Density, kg/m 3 

a Turbulence modeling constant 

t ij Reynolds stress tensor, Pa 

t w wall shear stress, Pa 

0 Momentum thickness, m 

English Symbols 

a 0 Turbulence modeling constant 

c Sound speed, m/s 

Turbulence modeling constant 
C p Pressure coefficient 

Cf Skin friction coefficient 

F\ Turbulence modeling constant 

F ‘2 Turbulence modeling constant 

k Turbulent kinetic energy 

M Mach number 

P Pressure, Pa 

V Turbulent model production term 

Re Reynolds number based on momentum thickness, m 

Rt Turbulent Reynolds number, m 

Sij Strain rate tensor, 1/s 

S t Stanton number, q w / pooUooC v [rT aw - T w ) 

T Static temperature, K 

u'v' Turbulent shear, m 2 / s 2 

u',v' Fluctuating velocity components, m/s 

x , y Cartesian body axes, m 

Subscripts 

aw Adiabatic wall 

eq Equilibrium condition 

k Kinetic energy equation constant 

o Total conditions 

w Wall conditions 

t Turbulent conditions 

oo Freestream conditions 

uj Specific dissipation equation constant 

Superscripts 

D Diviatoric part of a tensor 


I. Introduction 

Computational predictions of the aerodynamic and aerothermodynamic environments a space vehicle en- 
counters are essential for the design and operation of the vehicle. As computer speeds continue to increase, 
Computational Fluid Dynamics (CFD) can be used to not only predict the flowfield around simple config- 
urations, but also complex vehicles. The advances in computer speeds have given rise to advances in the 
modeling of turbulent flows using Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES). 
For complex configurations, DNS or even LES are still impractical because of the grid spacing requirements 
and the need for time-accurate solutions. Although DNS and LES employ a more physics based representa- 
tion of the fluid dynamics occurring, the Reynolds Averaged Navier Stokes (RANS) Equations solved along 
with a turbulence model are still a valuable tool for aerodynamic analysis. 

One difficulty with current one and two-equation turbulence models is the inability to account directly 
for non-equilibrium effects such as those encountered in large pressure gradients involving separation and 
shockwaves. Current turbulence models such as Spalart’s one-equation model, 1 the classic k — e 2,3 and 
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Wilcox’s k — uj 3 4 two-equation models have been designed and tuned to accurately predict equilibrium flows 
such as zero-pressure gradient boundary- layers and free shear layers. Application to more complex flows 
can be problematical at best. Although there have been many attempts to modify or correct basic one- and 
two-equation models, most of these attempts have been only marginally successful in predicting complex 
flows. 

More complex models such as Reynolds stress models have been investigated extensively, primarily for 
relatively simple flows but also for complex flows. In most cases these models give somewhat better pre- 
dictions than the simpler one and two equation models, but for complex flows they do not perform much 
better than the simpler models. Rizzetta 5 evaluated three popular explicit algebraic Reynolds-stress (EARS) 
models for separated supersonic flows, a shock impinging on turbulent boundary layer and a compression 
ramp. He compared these results to the k-e model. The results consistently showed that the EARS models 
offered little improvement over standard RANS approaches using the Boussinesq approximation. Viti et al. 6 
showed similar results for the Reynolds stress-transport turbulence models. 

One theoretical advantage of Reynolds stress models is that they directly account for non-equilibrium 
effects in the sense that the Reynolds stresses do not respond instantaneously to changes to the strain rate 
but more realistically lag them in time and/or space. Unfortunately, the Reynolds stress models are usually 
considerably more complicated and numerically stiff than the one- and two- equation models, and this has 
prevented their wide application for complex flows. 

In this paper we introduce a new turbulence model designed to account for non-equilibrium effects without 
invoking the full formalism of the Reynolds stress models. The idea is an extension of Olsen and Coakley’s 
lag model, ' which took a baseline two-equation model and coupled it with a third (lag) equation to model 
the non-equilibrium effects for the eddy viscosity. The third equation was designed to predict the equilibrium 
eddy viscosity in equilibrium flows. The proposed lagRST model uses the same lag equation idea, but directly 
lags the Reynolds stress tensor components instead of the eddy viscosity. The Boussinesq approximation 
sets the equilibrium values of each Reynolds stress tensor component that the lag equations will drive the 
value towards. 

We show results obtained with two Reynolds stress lag models based on the Wilcox k -u> model 4 and 
Menter’s SST model, 8,9 respectively. Applications to four flows are given including an incompressible flat 
plate flow, an incompressible mixing layer, 10 an essentially incompressible adverse pressure gradient flow 
with separation, 11 and a transonic bump flow 12 with a shock wave and separation. In addition, aerodynamic 
comparisons will be made to Apollo wind tunnel data. 13 We also are looking for a compressible shock wave 
turbulent boundary layer interaction case to assess. Results using the new model are compared with results 
obtained with Spalart’s model 1 and Menter’s k-u> SST model. Results obtained with the new model show 
encouraging improvements over results obtained with the other models. 


I. A. Non-Equilibrium Turbulence 

Several reviews have been done 14-16 assessing the accuracy of one- and two-equation turbulence models in 
the prediction of SWTBLI. The reviews show that one- and two-equation models poorly predict the starting 
position of the adverse pressure gradient flow, over predict the surface pressure within separation, and over 
predict skin friction and heat transfer at the re-attachment point. It is proposed that a main cause of 
this inaccuracy of the predictions in the separated regions is the effect of flow history. In other words, 
the Boussinesq approximation relates mean strain rates directly to the Reynolds stress tensor (RST). Thus 
any change in the strain rate is instantly accounted for in the RST. In reality, there is a time lag between 
changes in the strain rate and changes in the RST when the turbulence is not in equilibrium. Turbulent 
flows deviate from equilibrium conditions when large adverse pressure gradients are present. The standard 
one- and two-equation turbulence models have been designed and tuned to accurately predict equilibrium 
flows, such as zero-pressure gradient boundary layers and free shear layers. 

Clauser 1 ' defined a boundary layer to be in equilibrium if the pressure gradient parameter 


P 


S* dP 0 
t w dx 


(1) 


is a constant. When boundary layers are in equilibrium, they are self-preserving, i.e. when properly scaled, 
the wake portions of the velocity profile for two different boundary layers with the same value of f3 are 
identical, even if their Reynolds numbers are different. Equilibrium boundary layers have the characteristic 
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that the turbulent time scales are much smaller than the mean flow time scales so that the turbulence can 
react quickly to these slowly changing mean flows. 

I.B. lagRST Model 

Modeling nonequilibrium turbulent effects have been attempted by many researchers. An overview of many 
of these models can be found in Lillard. 18 One of these models, the Lag model of Olsen and Coakley' and 
Olsen et al., 19 has provided promising results for nonequilibrium flows. It is formed by taking a baseline 
two-equation model and coupling it with a third equation (lag equation) to model the nonequilibrium effects 
in the eddy viscosity. This class of models essentially introduces a lag into the response of the eddy viscosity 
to rapid changes in the mean flowfield, which is the character of the response seen experimentally. By lagging 
the third equation it provides the turbulence model an extra degree of freedom without affecting accurate 
equilibrium flow predictions. 

The lagRST model was designed to build upon the success of the Lag model by using the same base- 
line formulation, but directly solving for the Reynolds stresses, relaxing each stress independently along a 
streamline. The equilibrium state is defined by transport equations for the turbulent kinetic energy, k, and 
the specific dissipation rate, w. These are shown in equations 3 and 4, respectively. The lag equation for 
turbulent viscosity is then replaced with a similar relation for the Reynolds stress tensor, t, ,■ , defined as 

Tij = -pu'iV'y ( 2 ) 

Equation 5 details the lag equation for T,y. A separate equation is added for each Reynolds stress (this 
becomes six equations since the tensor is symmetric). 


where: 


dpk d 
dt + dx k 
dpu d 
dt dxk 


^ pu k k - [p + cr k pt] 
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with parameters 


( 3 ) 

( 4 ) 

( 5 ) 


«o = 0.35 $ = 0.075 

a = 5/9 p* = 0.9 


(jfc = 1.5 

cr e = 0.5. 


The equilibrium values of the Reynolds stresses are determined by using the Boussinesq approximation and 
solving for Tjj eq . 


T ijen ~ P 


2 (2 

— kSij v t eq ( 2sij -s kk 5ij 


( 6 ) 


I.B.l. Lagging the Reynolds Stress Tensor using Menter’s SST 

A version of the LagRST model using Menter’s SST model 9 to provide the equilibrium values of k and u 
was also implemented. Equations 7 and 8 provide the transport equations for the equilibrium quantities. 
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( 7 ) 
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(e + w)^- 
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0~ 


„ , 1 dfc <9w 
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Lv l/tl/j Ky tAs j 


(8) 


The modeling constants and blending functions used were not modified and therefore can be found in Menter. 9 
The turbulent eddy viscosity was calculated the standard way, Where v t is defined as 


aik 

teq max (aiw; 0 F 2 ) 

This was then substituted into equation 6 to obtain Tij eq . 


(9) 


I.C. Summary of Proposed Models 

Two Reynolds stress lag models were implemented into OVERFLOW v2.0aa and v2.2c. Table 2 details the 
name of the model and how the equilibrium value are obtained. The lagRST model uses the k — to model 
to define equilibrium values while the lagRSTSST model uses the SST model. Both models use the same 
transport equation for T,y, which was already shown in equation 5. The lagRST model per time step was 


Table 1. Lag variables 


Model (Designator) 

lag variable 

equilibrium variable 

LagRST (905) 

Tij = ~pK v j 

T ijeq = ~P {skeqd-ij " — 3 S kk Sij )) 

z/ t = — , defined from k — u 

‘-eq U e « ’ 

lagRSTSST (909) 

nj = -puy. 

v t = a ik eq defined from SST 

z *q max(aiuj eq ;ilF2) ’ 


35% slower than the SA model and 20% slower than the SST. 

I.C.l. Modeling Constants 

During the course of this study, it was discovered that an exhaustive study of the effect of a k had not been 
done in the original Lag model formulation. Since the formulation of the original Lag model, three values 
had been proposed for the value of a kl 0.5, 0.8 and 1.5. Although previously thought to only be a cosmetic 
change, the value of a k makes a substantial difference near the edge of the boundary layer. In addition, a 0 
will be assessed for values of 0.35 and 0.2. The value of 0.35 is the standard defined by Olsen and Coakley 7 
and Olsen et al. 19 

I.D. Numerical Method 

The lagRST and lagRSTSST models were implemented into version 2.0aa and 2.2c of the OVERFLOW 
code. OVERFLOW 20,21 is a structured (overset) grid, Navier-Stokes flow solver. It uses a finite-difference 
formulation, with flow quantities stored at the grid nodes. OVERFLOW has central- and Roe upwind- 
difference options, and uses a diagonalized, implicit approximate factorization scheme for time advancement. 
Local timestepping, multigrid techniques and grid sequencing are all used to accelerate convergence to a 
steady state. The standard turbulence models used by the OVERFLOW code are the one-equation SA 
model and the two-equation SST model. Other models available are the Baldwin-Lomax, k-e and the k -u. 
SA and SST are used for general aerodynamic applications within the OVERFLOW community. The full 
Navier Stokes equations were solved for all solutions. 
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I.E. Results with Similar Models 


Hamlington and Dahm 22,23 have derived a model similar to our proposed lagRST model in which they use 
the Reynolds anisotropy tensor (similar to the version discussed above). However, instead of numerically 
solving the lag equation, they solve the lag equation along a streamline, reducing it to an ordinary differential 
equation. They assessed several homogenous flows and have done work on flat plates and the interaction of 
an impinging oblique shock wave with a turbulent boundary layer. 

Churchfielcl and Blaisdell 24 used a similar definition of the lagRST model 25 discussed herein, but used 
the 2006 k — u> formulation to define vteq- Churchfield and Blaisdell 26,2 ' then went on to use a previous 
version of the lagRST formulation incorporated in the OVERFLOW code provided by Lillard 25 to assess 
wingtip vortex flow. 


II. Results 

Four test cases were chosen, to cover a broad range of flow physics while keeping the modeling complica- 
tions to a minimal level. These cases were chose in order to execute the following: 

1. Verify the implementation of the lagRST and lagRSTSST models. 

2. Confirm the value of the lag constant, a Q , previously derived by Olsen and Coakley 28 as 0.35. 

3. Assess the sensitivity of the models to ctk and confirm the value obtained by Olsen et al. 29 

4. Assess the performance of the lagRST and lagRSTSST models in predicting basic aerodynamic flows, 
such as turbulent boundary layer flow, mixing layers, and simple separated flows. 

In addition, several other turbulence models were run for comparison purposes including the Spalart- 
Allmaras, k -u) 1998, SST, and Lag models. Although the k — u model is not a standard model anymore due 
to known freestream dependencies, it will be assessed here for comparison because it is the base model for 
the Lag model and the majority of the lagRST models. Table 2 lists the set of turbulence models that are 
referred to by number instead of name in some of the upcoming figures. The associated number is also the 
turbulence model designator within the OVERFLOW code. 

Table 2. Table of turbulence models and their associated reference numbers in the OVERFLOW code 


Model Number 

Turbulence Model 

102 

SA 

202 

k -ijj 

203 

SST 

304 

Lag 

903 

Lag model with Reynolds tress tensor implementation 

905 

lagRST - Tij using k -w for equilibrium values 

908 

SST model with Reynolds tress tensor implementation 

909 

lagRSTSST - using SST for equilibrium values 


The four test cases chosen were the incompressible flat plate of Osterlund 30 and Osterlund et al., 31 
adverse pressure gradient boundary layer of Driver, 11 incompressible mixing layer of Bell and Mehta, 10 and 
a transonic bump flow by Bachalo and Johnson. 12 The comparisons shown in the following sections are part 
of a bigger study done by Lillard. 18 In particular, many results showing verification of code implementation, 
detailed grid convergence for multiple models, and velocity and Reynolds stress profile comparisons will not 
be shown here. 

II. A. Incompressible Flat Plate 

Osterlund 30 and Osterlund et al. 31 obtained zero pressure gradient turbulent boundary layer data in the 
Minimum Turbulence Level or Marten Theodore Landahl (MTL) wind-tunnel at the Department of Mechan- 
ics, Royal Institute of Technology (Swedish: Kungliga Tekniska Hogskolan, abbreviated KTH). The data was 
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taken on a 7 meter flat plate, with instrumentation locations ranging from 1.5 to 5.5 meters. Reg ranged 
from 2500 to 27000 with freestream velocities ranging from 10 to 50 m/s. Fluctuating velocity components 
and the fluctuating wall-shear stress were measured. Velocity data was taken with hot-wire anemometry 
probes mounted on a traversing system that protruded from the plate Skin friction was obtain by using both 
oil-film interferometry and using the mean velocity measurements and fitting them to the law of the wall 
with a least squares fit. This test case was first chosen to verify the following: 

1. The lagRST model’s ability to predict the near wall portion of the turbulent boundary layer (law of 
the wall). 

2. The model’s ability to predict the wake portion of the turbulent boundary layer. 

To compare velocity profiles, all the models were independently shifted so they matched the experimental 
value of 6 at x = 3.5 meters. This was done by changing the upstream running length of the grid (for each 
turbulence model) so that at x = 3.5 meters the predictions gave the experimental value of 0. This was 
necessary to remove transition effects from the model comparisons. The velocity profiles were then compared 
at x = 5.5 meters in the shifted coordinate space. 

II.A.l. Grid Convergence 

The grids generated were cartesian using geometric spacing in the off-body direction. The nominal grid 
chosen (referred to as the medium grid) had 92 streamwise points and 129 points in the boundary layer. 
The grid was extended well beyond the boundary layer edge to remove any boundary consdition effects. The 
wall spacing chosen gave a y + less than 0.15 across the plate. To create the coarse and fine grids, both grid 
dimensions were scaled by \/2 to increase the total grid dimension by a factor of 2. The wall spacing was 
then changed proportionally. This gave grids with dimensions of 65 by 91 and 129 by 182. 

Figure 1 shows the worst case for grid convergence for the lagRST model, with ag = 0.5. A non-physical 
inflection point near the edge of the boundary layer develops for this model. The inflection point begins 
around 2 = 0.04m. Not only is there a large deviation at the edge, but the non-physical behavior gets 
worse as the grid density progresses from coarse to fine. This behavior nearly goes away with ag = 0.8 and 
completely disappears with a g = 1.5. Table II.A.l presents the skin friction values near the end of the flat 
plate, Reg = 25,000, for the three variations in ag. The % difference listed is the difference for either the 
coarse or medium grid from the fine grid solution. ag = 0.5 clearly has the largest % differences, indicating 
a lack of grid convergence to go along with the non-physical velocity profile. 


Grid 

cr k = 1.5 

OO 

o 

II 

b 

ag = 0.5 

coarse 

0.002303 (0.35%) 

0.002263 (1.57%) 

0.02208 (2.75%) 

medium 

0.002300 (0.22%) 

0.002246 (0.81%) 

0.002176 (1.25%) 

fine 

0.02295 

0.002228 

0.002149 


Table 3. Skin friction values at Reg = 25000 for the lagRST model with a Q — 0.35 and varying with the % difference 
between the coarse or medium grid and the fine grid listed in parenthesis. 


II. A. 2. lagRST Constant Selection 

This section assess the six combinations of a 0 and ag and removes any pair that generates non-physical 
behavior. Figure 2 shows the lagRST model with varying coefficients compared to Osterlund’s velocity 
profile data. As it was seen before with previous plots, low values of a 0 and ag cause an inflection in 
the boundary layer profile that is not physical. The profiles that retain a shape without an inflection are 
a 0 = 0.35 along with ag = 0.8 and both solutions with ag = 1.5. Although clearly not an exact match the 
the experimental data, the pair of a 0 = 0.35 and <Jk = 0.8 provides a good representation without having 
the inflection point the lower values values are plagued with. 

Based on the previous figure, several combinations of the constants can be discarded. The inflection in 
the velocity profiles removes ag values of 0.5 and low a 0 values with ag = 0.8. The skin friction comparisons 
reinforce the removal of a Q = 0.2 with erg = 0.5 due to how much lower they are than the correlations and 
data. That leaves three sets of coefficients to assess, a 0 = 0.35 with ag = 0.8 and a 0 = 0.2 or a 0 = 0.35 with 
ag = 1.5. 
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U(m/s) 


Figure 1. Velocity profiles for lagRST with a & = 0.5 showing the effects of varying grid dimensions 



U(m/s) 


Figure 2. Velocity profiles for the lagRST compared to Osterlund’s data. 
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II. A. 3. Model Comparisons 


Figure 3 shows the lagRST model with the previously selected combinations of a 0 and <Tk as well as the 
lagRSTSST model compared to the standard SST and SA models. Unfortunately, none of the models match 
the data across the entire profile. Although none of the models are perfect, it is worth noting that the 
general error is all on the order of 5 — 10%. In addition, the inflection near the edge in the lagRST models 
is avoided based on the selection of constants. For the lagRST models, the two highest coefficients have the 
most deviation from the experimental data. It can also be said that the edge of the boundary layer the SA 
model under predicts the velocity magnitude more than the majority of the models. 



U(m/s) 


Figure 3. Velocity profiles comparing the Lag, lagRST, lagRSTSST, SST, and SA models 


Figure 4 shows the same results in law of the wall space with Coles law of the wake added for comparison. 
The SA model predicts the transition from the laminar sublayer to the law of the wall region much more 
accurately than all other models. However, the slope in the law of the wall region is slightly in error compared 
to the other models. On this scale, the SST models gives a very close result to the Lag and lagRST models 
with constants a o and a j. equal to 0.35 and 0.8, respectively. None of the models predict the extent and 
magnitude of the wake region properly. 

Skin friction predictions for all the models were within 2% of each other and are therefore not shown 
here. 


II. B. Mixing Layer 

Bell and Mehta 9 10 obtained experimental data for a two-stream mixing layer with a velocity ratio of 0.6 in 
a Mixing Layer Wind Tunnel located in the Fluid Mechanics Laboratory at NASA Ames Research Center. 
Measurements were taken with a rotatable cross-wire probe held with a three-dimensional traverse. They 
obtained data for both tripped and untripped boundary layers. For this comparison, only the tripped 
turbulent data will be used. They defined the shear-layer thickness, S, using a least-squares fit of the mean 
data to the error function profile shape: 


where r) is defined as 



(10) 
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Figure 4. Velocity profile in law of the wall space comparing the Lag, lagRST, lagRSTSST, SST, and SA models with 
the addition of Cole’s law of the wake. 


and U* , the velocity parameter, is defined as 


V = 


Y-Y 0 

S 


( 11 ) 


u* 


u-u 2 

t/l-t/2 


(12) 


and Y a is the centerline of the mixing layer and is defined from the error function fit. The average RMS error 
for the least squares fit had less than a 2% deviation between the functional form and the data. To compare 
to this data, the velocity profiles were extracted and then processed through a least squares fit to calculate 
<5 and Y a . The upper velocity stream, U upper , was set to 15 m/s while the lower velocity stream, Ui owe r, 
was set to 9.0 m/s. The solution was initialized for all positive z values to U upper and negative z values to 
Uiower- The inflow condition was set to the same velocity condition as the solution initialization and used an 
inflow/outflow condition. The upper and lower conditions were set to U upper and Ui ower , respectively and 
used the same inflow/outflow condition as the initial condition. The downstream boundary condition was 
set to outflow (extrapolation). 


II.B.l. Grid Convergence 

Three grids were generated in order to study grid convergence for this test case. The finest grid ran was 
129 (axial) by 205. The grid was then decreased to 65 by 103 and further to 33 by 52. Figures 5 shows 
grid convergence for the lagRST and lagRSTSST models. For the lagRSTSST model, the medium grid (65 
by 103) gets within 2% of the fine grid answer. For the lagRST again the medium grid gives excellent grid 
convergence and is within 1%. 


II. B. 2. Model Comparisons 

Figure 6 shows a comparison of the SST, lag, lagRST, and lagRSTSST models. The lagRSTSST compares 
very well with the standard two-equation SST model and the experimental data. The lagRST and Lag model 
comparisons are very poor, even with large variations in <jfc. To date, there has been no explanation for this. 
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x(m) 


Figure 5. Grid convergence for the lagRST and lagRSTSST models compared to the mixing layer growth data for the 
fine (129 by 205), medium (65 by 103), and coarse (33 by 52) grids. 


For the lagRSTSST model, a lower value of the lag constant, a 0 = 0.2, was assessed for this case. It had 
little effect, but did increase the slope of the growth rate by 1.5%. 

II.C. Driver’s CSO Flowfield 

Driver 11 executed an incompressible test of a turbulent separated boundary layer in a low speed open circuit 
wind tunnel. The test geometry was a circular cylinder aligned with the freestream and mounted in the 
center of the tunnel. It was surrounded by flexible tunnel walls that could be diverged in order to create 
an adverse pressure gradient. Skin friction measurements were taken with an oil-flow laser interferometer. 
Profile data was taken with a three-component laser Doppler velocimeter (LDV) system. 

The computational geometry was modeled with an axisymmetric cylinder to define the surface and an 
external streamline taken from the experimental data to define the opposite boundary. The outer streamline 
was treated as an inviscid wall. The viscous wall is no-slip and adiabatic. The upstream boundary condition 
was defined as constant total pressure and temperature, with static pressure allowed to vary and velocity 
direction aligned with the cylinder axis. The downstream static pressure was adjusted at the exit to match 
the experimental static pressure upstream of the interaction region ( x = —0.438m). This was an iterative 
process that was done for each simulation. In addition, the upstream length was also taken from Olsen and 
Coakley 7 to match the experimental boundary layer thickness. 

The calculations were made on the same grids used in the original analysis by Olsen and Coakley.' The 
fine grid had 200 (axial) by 160 (radial) points. The coarse grid had a factor of 2 reduction in grid points 
(100 by 80) and the medium grid had 150 by 120. Olsen used the grids referred to as fine and coarse in his 
study. Figure 7 shows the computational domain for this test case. 

II.C.l. Grid Convergence 

Figure ?? shows the skin friction distribution for the lagRSTSST model. The model has excellent convergence 
behavior. There are mild differences well upstream, but once the adverse pressure gradient area starts all 
models are in agreement. The pressure comparisons were nearly identical between the grids, thus they are 
not pictured here. 

The medium and fine grids for all cases matched the measured upstream boundary layer thickness to 
within 10%. 
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Figure 6. Mixing layer growth rate comparisons for the SST, lag, lagRST, and lagRSTSST models. 



Figure 7. Computational domain, coarse grid pictured. 
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Figure 8. 


Skin friction comparisons showing grid convergence for the lagRST model with standard coefficients. 


Although now shown here, the lagRST model gave comparable results. The medium and fine grids 
for all cases matched the measured upstream boundary layer thickness to within ten percent. It was also 
confirmed that the grid system had acceptable convergence performance when looking at the Reynolds 
stresses (although not shown here). 

II. C. 2. Model Comparisons 

In previous sections, the lagRST coefficients have been reduced to three sets of choices. For ease of reference, 
table II. C. 2 shows the new naming convention for the lagRST based models for the three different sets of 
coefficients. 


Table 4. Model naming convention for lagRST models 


lagRST model name 


a a 

lagRST-1 

0.8 

0.35 

lagRST-2 

1.5 

0.2 

lagRST-3 

1.5 

0.35 


Figure 9 shows the pressure distribution versus axial distance for the SST, Lag, lagRST, and lagRSTSST 
models. For the lagRST models, the lagRST-2 model matches the pressure more accurately. However, all 
the lagRST models are an improvement over the baseline k — u> model. The Lag and lagRST-2 provide the 
best prediction for these models. The SST-based models provide the best overall prediction of pressure. 
The lagRSTSST model predicts the pressure to have a lower plateau in the interaction region than the SST 
model. Overall, adding the lag equation improves the predictions for the baseline models. The lagRSTSST 
model was assessed with a lag constant of a Q = 0.2 for this case. It had less than a 1% change on the size 
and location of the separation zone vs the lagRSTSST model with the standard coefficient. 

Figures 10 and 11 detail the skin friction distribution on the cylinder. There is insufficient data to 
exactly characterize the separation zone length using skin friction, so the experimental separation length was 
determined with oil flow interferometry. The separation zone ranged from approximately from x = 30 mm 
to x = 220 mm. This was within ±25 mm of the mean flow data. Table II. C. 2 lists the separation and 
reattachment points as well as separation length for all models. 
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Figure 9. Pressure versus streamwise locations for the Lag, lagRST, and lagRSTSST models compared to k — uj and 
SST. 



x(m) 


Figure 10. Skin friction versus streamwise locations for the Lag, lagRST, and lagRSTSST models compared to k — lj 
and SST, with the separation noted by arrows. 
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Figure 11. Skin friction versus streamwise locations for the Lag, lagRST, and lagRSTSST models compared to k — lj 
and SST focusing on the separation region. 


All the k — u) based RST models provide a more accurate separation zone length prediction than the 
SST-based models. By far the best comparison comes form the lagRST model with ao = 0.35 with either 
value of <Tk chosen (lagRST-1 or lagRST-3). It is interesting to note that the baseline k — u> model is the only 
one which predicts a smaller separation zone length than the experiment. The lagRSTSST model provides 
the worst separation zone length prediction of any model. The majority of the models do a better job 
characterizing the reattachment point than the separation point. None of the models provide a prediction 
of the separation point that is within 12% of the experimental results, whereas two of the models provide 
reattachment predictions within 12% of the experimental results. 

Table 5. Table of data indicating location and extent of the separation zone with percent error compared to the tunnel 
data. The percent error uses the separation zone length in the denominator. 


Model 

Separation Point 

Reattachment 

Zone length 

Data 

0.03m ± 0.025 

0.22m ± 0.025 

0.19m ± 0.025 

k-w 

0.00708 (-12.1%) 

0.131 (-46.6%) 

0.124 (-34.5%) 

SST 

-0.015 (-23.6%) 

0.258 (20.0%) 

0.273 (43.6%) 

Lag 

-0.0269 (-29.9%) 

0.226 (3.16%) 

0.253 (33.1%) 

lagRST-1 

-0.011 (-21.6%) 

0.199 (-11.2%) 

0.210 (10.4%) 

lagRST-2 

-0.00955 (-20.8%) 

0.233 (6.67%) 

0.242 (27.5%) 

lagRST-3 

0.00328 (-14.1%) 

0.197 (-12.2%) 

0.194 (1.8%) 

lagRSTSST 

-0.0168 (-24.6%) 

0.259 (20.4%) 

0.276 (45.0%) 


II. D. Transonic Bump Flow 

Bachalo and Johnson 12 tested a transonic interaction on an axisymmetric cylinder with a circular arc bump. 
This test was executed in the NASA Ames Research Center 2x2 foot Transonic Wind Tunnel. Data was 
taken with a two velocity-component laser velocimeter. The upstream boundary conditions were constant 
total pressure and temperature, with static pressure allowed to vary and velocity direction aligned with the 
cylinder axis. The outer edge of the flowfield was extended approximately 10 bump cord lengths away from 
the cylinder surface and utilized a characteristic boundary condition. The viscous wall was no-slip adiabatic 
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with the downstream static pressure held at pi n f. To maintain consistency with Olsen and Coakley, 7 the 
upstream length of the cylinder was maintained at the value they used, which provided a boundary layer 
thickness of approximately 1 cm at the location of the bump, which was the same approximate thickness 
measured in the test. 



x/c 


Figure 12. Pressure plots for the lagFtST-l model showing grid convergence. 


II.D.l. Grid Convergence 

Three grids were chosen for this case. A coarse grid with dimensions of 96 (axial) by 101 (off-body), medium 
grid with dimensions 192 by 151, and a fine grid with dimensions 390 by 201. Wall spacings that translated 
to y + values less than 0.2, 0.1, and 0.05 were used for the coarse, medium, and fine grids, respectively. 

Figure 12 shows the pressure comparisons for the lagRST-1 model for the various grid resolution cases. 
The medium and fine grid give nearly identical results. The other models (lagRST-2, lagRST-3, and la- 
gRSTSST), all had similar convergence behavior, if not better. The medium grid provides adequate grid 
convergence for this test case. 

II. D. 2. Model Comparisons 

Figure 13 shows the pressure comparisons for the lagRST based models as well as the k -u>, SST, and Lag 
models focused in on the the interaction region. The prediction of the shock location is of vital importance, 
but is not consistently predicted by the models. The k -oj predicts a very late shock location and the 
lagRSTSST predicts it very early. These two models bracket the other models. The SST and the lagRST-2 
models give nearly identical answers. They predict a slightly upstream location for the shock. The Lag and 
lagRST-1 also give nearly identical answers to each other, and provide a very accurate prediction for the 
shock location. The lagRST-3 also gives a very accurate shock location. None of the models predict the 
right pressure behavior in the separation zone (a; = 0.71m — 1.1m). The lagRSTSST model has the closest 
comparison. However, none of the models predict the steep drop-off in pressure just after re-attachment that 
is present in the experimental data. 

Besides shock location and pressure comparisons, perhaps the more important comparison to make is of 
the separation zone length. Table II. D. 2 shows the observed separation and attachment points as well as 
separation zone lengths for each turbulence model. The table indicates the k-w model provides an accurate 
prediction of the separation zone, but this only refers to the extent of the separation. Previous figures have 


16 of 37 


American Institute of Aeronautics and Astronautics 



Figure 13. C v versus axial location for all turbulence models compared to the experimental data focusing on the 
separation zone. 


shown that the pressure prediction in the separation zone from this model is under predicted. The Lag and 
lagRST-1 models provide the best characterization of the separation zone length, with the worst comparisons 
coming from the SST and the lagRSTSST models. The lagRST-3 and lagRST-2 models (which both use a 
larger value of a*,), also have a large error in the prediction. The lagRSTSST model was assessed with a lag 
constant value of a Q = 0.2. This had no effect on separation and less than 0.1% change on the separation 
zone length and position. 

Table 6. Table of data indicating location and extent of the separation zone with percent error compared to the tunnel 
data. The percent error uses the separation zone length as the denominator. 


Model 

Separation Point 

Reattachment 

Zone length 

Data 

0.71 

1.1 

0.39 

k -u 

0.701 (-3.03%) 

1.05 (-13.8%) 

0.346 (-10.8%) 

SST 

0.645 (-17.4%) 

1.18 (21.6%) 

0.539 (39.0%) 

Lag 

0.666 (-11.9%) 

1.12 (4.64%) 

0.452 (16.6%) 

lagRST-1 

0.668 (-11.4%) 

1.12 (4.14%) 

0.210 (10.4%) 

lagRST-2 

0.646 (-17.1%) 

1.19 (21.9%) 

0.539 (39.1%) 

lagRST-3 

0.659 (-13.7%) 

1.15 (12.5%) 

0.5 (26.3%) 

lagRSTSST 

0.632 (-20.8%) 

1.19 (23.3%) 

0.556 (44.1%) 


II. E. Shock Impinging on a Turbulent Boundary Layer 

Schulein 32 ’ 33 conducted a = 5) wind tunnel test in the DLR Gottingen Ludwieg Tube facility, using air 
as the test gas. The nominal freestream conditions for this test are listed in table 7. The measured values 
included wall pressure, skin friction (by an oil-flow technique), and wall heat transfer (infrared camera 
measurements and thin film gages). The provided experimental uncertainties were approximately 2% for 
pressure, 10% for skin friction, and 20% for heat flux. The upstream boundary layer was fully developed 
and an assessment by Brown 34 showed that it adhered to standard turbulent correlations. 

Schulein tested four configurations, starting from a 2-D nominal zero-pressure gradient flat plate boundary 
layer and then adding a shock generator with 3 different angles, 6°, 10°, and 14°. Brown 34 did an exhaustive 
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Table 7. Nominal test conditions for Schulein’s experiment. 


Condition 

Value 

Pt 

T 0 

T w 

Re/m 

Coo 

Rq,oo 

2.12 MPa 
410K 
300K 
36zl0 6 
830 f 
0.41# 


study of all three angles and provided recommendations for which cases to run. Other comparisons to this 
dataset are available in the literature, including Fedorova et al. 35 and Steelant. 36 For this study, the 2-D 
boundary layer was assessed to ensure the quality of the numerical upstream profile before the interaction. 
The only shock generator angle assessed was the 14° case. This case has the largest separation zone tested 
by Schulein. In addition, it has been chosen by NASA’s Fundamental Aeronautics Program as a major test 
case to use for model comparison. The lower shock generator angles have attached or incipient separation, 
thus do not assess the separation behavior of the model. 

An issue that plagues most cases that have this configuration is the shock generator is not relatively 
“long” compared to the nominal boundary layer thickness. This allows the expansion waves that emanate 
from the tail of the shock generator to interfere with the separation region, making the interaction more 
complex and not allowing for the various topological features of the SWTBLI flowfield to develop. Schulein 
was able to have a ~ 1000, which is a large enough ratio to provide ample space for the interactions to 
set up. 



Figure 14. Coarse grid for the 14° SWTBLI case (blue line indicates the shock generator). 


II.E.l. Grid System 

The initial grid system was identical to the set run by Brown. 34 Table 8 gives the dimensions and y + values 
for each grid system provided by Brown. Figure 14 shows the grid for the 14° SWTBLI, where flow would 
be from left to right. The blue line indicates the viscous wall of the shock wave generator. The adjacent 
boundary conditions to the shock generator were set to inviscid walls, and the opposite wall was set to a 
viscous, isothermal wall. 


Table 8. Grid dimensions and y + values for the nominal set of grids. 


Grid System 

streamwise 

off-body 

y + 

coarse 

529 

129 

0.1 

medium 

1057 

257 

0.05 

fine 

2113 

513 

0.025 
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II. E. 2. Flat Plate Boundary Layer Results 


This case followed the standard method used in all other cases to setup the incoming turbulent boundary 
layer profile. This profile will be generated on the actual geometry by transitioning the simulation in order to 
match the momentum thickness calculated from the experimental velocity profiles. This is not the method 
Brown 34 used, as he numerically set the transition point based on the experimentally observed location. This 
method did not work in OVERFLOW, as there is no way to control the transition zone length. In particular, 
each turbulence model behaves differently in this transition zone, which made the results vary based on the 
model. Therefore, the previous method of setting the transition point to match the experimentally provided 
0 was used. This method provides the best estimate of the incoming turbulent boundary layer. 



x(m) 


Figure 15. Momentum thickness verses axial distance for the SST and lagRST-l models compared to Schulein’s 2-D 
boundary layer data. 


Figure 15 shows the matching process for the SST and lagRST-l models. The dashed lines are the fully 
turbulent cases while the solid lines have been tripped to match 9. Both cases were transitioned at different 
axial locations. The same process was done for all other turbulence models, but will not be shown here. This 
process shifts the virtual origin of the simulations so that the boundary layers between different simulations 
and the experimental data are in the same reference frame. Figure 16 shows the heat flux for the same profiles. 
The predictions go through the scatter in the data. This profile shows that once sufficiently downstream, 
the tripping process makes little difference on the value of heat flux (shown by the small differences between 
the dashed and solid lines of the same color) . The Reynolds number is sufficiently high that downstream of 
the interaction, the transitional effects on heat ffux are minimal. 

II.F.3. Grid Convergence 

The baseline set of grids run for this case were provided by Brown. 34 The dimensions were previously listed 
in table 8. All grid resolution cases were run fully turbulent to avoid sensitivities to transition location 
and the extent of the transition zone. Figure 17 shows the pressure distribution for all three grids for the 
SST model. The medium and fine grid show excellent grid convergence, even for the separation location. 
The lagRST model did not converge with the medium grid, but further testing showed the fine grid was 
converged. 
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Figure 16. Heat flux (W/m 2 ) verses axial distance for the sst and lagR.ST-1 models compared to Schulein’s 2-D boundary 
layer data. 
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Figure 17. Pressure ( Pa ) distribution for the SST model for the coarse, medium, and fine grid. 
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II. E. 4- 14° SWTBLI Data Comparisons 


Table II. E. 4 lists the separation and reattachment point along with the separation zone length for each model 
run. The SST model has by far the best characterization of the separation zone length and location. The 
percent errors between the simulations and the experiment for the separation zone length are less than 1%. 
The remaining models do a very poor job predicting the separation point, with the majority of the models 
having percent errors over 30%. The reattachment point is well characterized by all the models (numerical 
prediction of reattachment location is within 0.1% for all the simulations). The separation zone length 
thus has the same characteristics as the separation point uncertainties. The SST provides an extremely 
accurate answer and the remaining models have percent errors greater than or equal to 30%. The accuracy 
of the simulations at reattachment is an interesting result. Previous computations with the Lag model 
by Oliver et al. 37 39 on similar SWTBLI (Mach 3 compression ramps) had variations in the prediction of 
both the separation and re-attachment points. It is not clear why the shock impinging simulations predict 
the experimentally observed re-attachment point. The re-attachment appears to be driven by the external 
streamlines being processed through multiple shocks which re-align the inviscid flowfield with the viscous 
wall. However, this observation is largely geometric in nature, thus the exact physics of this phenomenon 
are unknown at this time. 

Table 9. Table of data indicating location and extent of the separation zone with percent error compared to the tunnel 
data. The percent error uses the separation zone length as the denominator. 


Model 

Separation Point 

Reattachment 

Zone length 

Data 

0.314m 

0.348m 

0.034m 

SA 

0.327 (39.2%) 

0.349 (4.2%) 

0.0221 (-35.0%) 

SST 

0.314 (0.6%) 

0.348 (0.7%) 

0.034 (0.1%) 

Lag 

0.300 (42.5%) 

0.349 (3.5%) 

0.0496 (46.0%) 

lagRST-1 

0.296 (52.9%) 

0.351 (7.7%) 

0.0546 (60.7%) 

lagRST-3 

0.301 (39.5%) 

0.348 (7.0%) 

0.0498 (46.6%) 

lagRSTSST 

0.304 (29.1%) 

0.348 (0.7%) 

0.0441 (29.8%) 

lagRSTSST low a 0 

0.295 (56.6%) 

0.35 (4.9%) 

0.0549 (61.5%) 



x(m) 


Figure 18. Pressure (Pa) distribution comparisons for all turbulence models on the fine grid. 
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Although the SST model matched the separation zone extent nearly perfectly, it does under predict 
the pressure plateau in the separation zone. Figure 18 compares the wall pressure distribution for all the 
models. The SA and SST models clearly miss the pressure plateau by around 30%. This is in contrast to the 
lagRSTSST and the lagRST models which match the pressure plateau very well. The lagRST-1, lagRST-3, 
and lagRSTSST with the low a Q over predict the pressure plateau, but by a small percentage. The under 
predicted pressure in the separation zone is contradictory to the observation of Wilcox 14 previously discussed. 
Wilcox found that the pressure in the separation zone was typically over predicted by turbulence models. 
The opposite trend in pressure comparisons is true for the pressure peak behind the reflected shock. The SA 
and SST models predict the peak very well, whereas all the lag-based models under predict the peak across 
the plateau region until the very end of the region. The SA and SST models predict the very flat pressure 
profile measured in the wind tunnel data. 



x(m) 


Figure 19. Skin friction distribution comparisons for all turbulence models on the fine grid. 


Figure 19 shows the skin friction distributions for all models run. Just as the comparisons to the undis- 
turbed boundary layer showed (figure ??), the skin friction is well under-predicted by all models. It is worth 
noting that Brown 34 presents the same under-prediction for this case. Brown used a second method to 
calculate the skin friction, which brought the skin friction results down by 20% - 30%. He analyzed the 
log-law region of pitot probe surveys, which he claims provides an equally accurate result in areas where the 
flow is attached. For this particular case, he believes that there may be instabilities in the oil flow due to the 
high shear rates in this area, which would lead credence to the delta between the two methods. The green 
dots show these results, which do have excellent agreement with the SST model. In this plateau region, the 
SST model has the largest deviation from the other models. The lag-based models all have the same shape, 
with the lagRSTSST having a slightly higher skin friction than the k — uj based models, or the lagRSTSST 
with the lag constant, a 0 , set to 0.2. 

Figure 20 shows the heat flux comparisons for all the models. All of the models over predict the heat 
flux by at least 25%, if not substantially more. The SST model has the worst behavior just downstream of 
the reflected shock, showing an extremely large overshoot in this region. The lag-based models have this 
overshoot, but the magnitude is much smaller. The SA model does not have this overshoot, but still over- 
predicts the heat flux. One important item to notice from figure 20 is that the lagRSTSST model largely 
removes the overshoot of the SST model. In fact, towards the end of the plateau region, the two models 
predict the same heat flux. So although lagging the Reynolds stress negatively affects the separation zone 
prediction, it makes a large improvement on the over prediction of heat flux. Note the heat flux measurements 
have an observed uncertainty of ~ 20%. 
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Figure 20. 


Heat flux ( W/m 2 ) distribution comparisons for all turbulence models on the fine grid. 


II. F. Orion Capsule 

As part of the CEV (now Orion) program, an experimental study 40,41 of the CEV geometry (shown in 
figure 21) was conducted at the Ames Unitary Plan Wind Tunnel (UPWT) 42 in the 11-foot Transonic Test 
Section and the 9X7 foot Supersonic Test Section. Test data were obtained for conditions at 0.3 < < 2.5, 

140° <a< 170°, and l.Ox 10 6 < Ren < 5.3xl0 6 . The angle of attack was focused around the trim condition. 
Three models were constructed during this test, one 7.66% scale and two 3.0% scale models. The larger 
model included pressure taps, unsteady pressure transducers, and an internal balance. One 3.0% model 
was designed exclusively for force and moment data (internal balance), and the other was designed solely 
for surface pressure data. The 3% models provided Reynolds number comparison, and allowed comparison 
between test run in two different facilities, which helped quantify uncertainty. The 7.66% scale model 
had both pressure taps and an internal balance. Boundary layer trips were applied to the model using 
CadCut trip dot tape. This study simulated the Ren = 5.3 x 10 6 case at Moo = 0.95. This condition 
was investigated only with the larger model in the 11-foot Transonic Test Section. Three different angles 
of attack were assessed, 142.4°, 154.4°, and 170.4°. Note that a = 180° would have the heatshield into the 
wind and a = 0° would be the cone into the wind. Through unpublished results and internal reports, the 
Orion Aerosciences team has shown that the choice of turbulence model makes a large difference on the wake 
behavior, which directly affects the Cd of the capsule. Therefore, this study will focus on comparisons with 
Cd and Cl- In addition, pitch plane contour images are used to show the wake size and characteristics. 

II.F.l. Grid System 

The grid system used was the same grid system employed by Olsen et al. 29 The capsule is modeled by a 
three zone grid system, shown in figure 22. The near body grid system is made up of an axi-symmetric grid 
(in y-z) that is rotated around the x-axis and two H-grids that are overset on the axis points to remove the 
singularities from the grid-system. The standard overset grid methodology within OVERFLOW 43 was used 
to generate the near body grids. The maximum y + of the grid systems were 1.5, 0.7, and 0.35 for the coarse, 
medium, and fine grids, respectively. The near field grid system extended out approximately from the 
capsule. The near body grids were then enclosed with cartesian box grids using the off-body grid generation 
capability internal to OVERFLOW 2.2c. 44,45 The wake box grid that enclosed the near body grids extended 
upstream of the capsule a little less than two capsule diameters and downstream over five capsule diameters. 
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Figure 21. Outer mold line geometry 


Heat-shield shoulder 
(A0 = 3 degrees) 


Nose shoulder 
(min_cell = 7) 



Figure 22. Near body grid system 
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This wake box grid was constructed of equispaced cells in all directions. The coarse, medium, and fine grids 
has A values of 10, 5, and 2.5 inches, respectively. The automatic techniques generated 30, 36, and 45 box 
grids for the coarse, medium and fine grids, respectively. An image of the medium grid system is shown 
in figure 23. The domain extended approximately 16 capsule diameters in each direction. The total grid 
dimensions are listed in table 10. 


Table 10. Number of grid points for the Orion grid system. 


Grid System 

Capsule Grids 

Box Grids 

Total 

coarse 

medium 

fine 

230,352 

1,788,565 

13,402,440 

1,014,386 

4,636,564 

32,323,258 

1,244,738 

6,425,129 

45,725,698 



Figure 23. Cartesian box grids enclosing the Orion near body grid system. 


II. G. Solution Procedure 

Simulations were initialized using multigrid and grid sequencing and completed with a time accurate sim- 
ulation using dual time integration to assess the unsteadiness of the flowfields. The time histories of the 
integrated forces were used to determine whether the flowfields were steady or unsteady. The integrated 
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forces also guided judgement about convergence and solution quality. Previous work by Olsen et al. 29 
showed that in general, supersonic cases were steady, and subsonic cases unsteady. This statement came 
from comparisons for a set of supersonic cases (Moo = 1-4) and a set of subsonic cases ( M oo = 0.95). Since 
this study will concentrate on the M a 0 = 0.95 cases, all cases will be simulated assuming unsteady flow 
conditions. 

II.G.l. Grid Convergence 

Due to the significant amount of computer time taken to run the cases, grid convergence was not assessed for 
every turbulence model. Grid convergence was assessed for the lagRST-1 model at all three chosen angles 
of attack, while the lagRSTSST model was only assessed at 154.4° and 170.4°. 

Table 11. Cl and Co for the a. — 170.4° set of cases. 


Model 

grid 

Co 

C L 

lagRST-1 

coarse 

1.10 

0.18 

lagRST-1 

medium 

1.08 

0.17 

lagRST-1 

fine 

1.09 

0.18 

lagRSTSST 

coarse 

1.22 

0.20 

lagRSTSST 

medium 

1.19 

0.19 

lagRSTSST 

fine 

1.19 

0.20 


Table 11 shows the Cl and Co for a = 170.4°. The results show for both the lagRST-1 and the 
lagRSTSST models that grid convergence is obtained with the medium grid. Even the coarse grid results are 
nearly identical to the medium grid. It is expected that this is due to both models producing steady solutions 
at this condition. The experimental results provided unsteady wake results, but none of the computations 
exhibited this behavior. This will be discussed in upcoming sections. The extent of the wake grid was also 
assessed by increasing its downstream distance by 20%. This had less than a 1% change on Co or Cl and 
did not affect the unsteadiness of the wake. 

Table 12. Cl and Cd for the a. = 154.4° set of cases. 


Model 

grid 

Co 

C L 

lagRST-1 

coarse 

1.09 

0.49 

lagRST-1 

medium 

1.06 

0.46 

lagRST-1 

fine 

1.04 

0.46 

lagRSTSST 

coarse 

1.06 

0.49 

lagRSTSST 

medium 

1.03 

0.46 

lagRSTSST 

fine 

1.01 

0.45 


Table 12 shows the Cl and Co for a = 154.4°. Again, the medium and fine grid results are nearly 
identical, but for this case, the coarse grid results have more variation. The lagRST-1 model did produce an 
unsteady result, whereas the lagRSTSST model produced a steady result. The medium grid provides a grid 
converged result. 


Table 13. Cl and Cd for the a. = 142.4° set of cases. 


Model 

grid 

Co 

C L 

lagRST-1 

coarse 

1.02 

0.65 

lagRST-1 

medium 

0.87 

0.59 

lagRST-1 

fine 

0.88 

0.60 


Table 13 shows the Cl and Co for a = 142.4°. Due to computing constraints, the lagRST-1 model was 
the only one used for this a. The coarse grid results for this case are very different than the medium and fine. 
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The medium and fine grid provide very similar results. This will be discussed more in upcoming sections, 
but the large deviation from the coarse grid is most likely due to the solution becoming more unsteady as the 
a gets lower. The numerical oscillations were larger at this condition, and the experimental results showed 
more unsteadiness. Based on this set of comparisons, the medium grid provides a grid converged solution 
and will be used for comparisons to the experimental data. 

II. G. 2. Flowfield Characteristics 

Results were obtained on the medium grid at the three previously discussed a conditions for a suite of 
turbulence models, ranging from baseline models to the lag-based models. The standard OVERFLOW 
versions of SA and SST were run along with the standard Lag model. In addition, the lagRST-1, lagRST- 
2, and lagRST-3 models were run along with the lagRSTSST model with the normal lag constant value 
(■ a 0 = 0.35) and a lower value of the lag constant, a 0 = 0.2. All results were averaged over a time interval 
containing a minimum of two wake oscillations. If an oscillation was not identified (either because the result 
was steady or the oscillation was not periodic), the time averaging interval was 3000 iterations. With the 
chosen time advancement parameters, a typical oscillation of the wake took at most 350 iterations. 



Iterations 


Figure 24. Time history of C D (plotted vs iteration) for the lagFtST-l model at all o’s. 


Figure 24 shows a typical time history of Cd for the lagRST-1 model at each a (the lagRST-3 model 
had similar behavior). The oscillation is periodic, indicating a consistent shedding behavior in the wake. 
The frequency of the oscillation increases going from a = 142.4° to a = 154.4°. The amplitude of the 
oscillation goes to zero as the a approaches 170.4°. Once at the highest a, the oscillation damps out and 
the Cd results become steady. Only a subset of the models had an oscillating Cd at the lower a’s, but 
every model produced a non-oscillating Cd profile at a = 170.4°. Figure 25 shows Mach contours for all 
three a conditions for the lagRST-1 model. In general, an oscillating Cd profile correlates to a wake that 
is shedding. This can be seen by comparing figures 25(a) and 25(b) with figure 25(c). At a = 170.4°, the 
Cd time history is steady, and the wake is symmetric and has no shedding. At the two lower a’s, the Mach 
contours show an asymmetric wake that is shedding, which correlates to previously seen oscillating Cd time 
history. This relationship will generally be the case for all solutions. If the wake is shedding, the Cd time 
history will be oscillating. There is however, one exception to this. Figure 26 compares the Lag model with 
the lagRST-1 model at a = 142.4°. The Lag model has a nearly steady Cd time history (variation in Cd 
of less than 0.5%), but figure 26(a) clearly shows the wake is oscillating. If you compare figure 26(a) with 
figure 26(b), it is apparent that the re-circulation zone in the wake of the capsule is closing for the Lag model 
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before the shedding is occurring. In contrast, the lagRST-1 simulation does not predict the re-circulation 
zone to close before the shedding begins. It stands to reason that if the re-circulation zone is closed, then 
the Cd time history will be nearly steady even if the wake oscillates downstream of the zone closure. 



(a) 142.4° 


(b) 154.4° 


(c) 170.4° 
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(d) Mach Contour Legend 


(e) Mach Contour Legend 


(f) Mach Contour Legend 


Figure 25. Mach contours in the pitch plane for the lagRST-1 model showing the wake at all three o’s. 


In contract to the lagRST-1 model, the lagRST-2 model shows an unsteady but non-periodic behavior 
in Cd- Figure 27 shows the Cd time history for over 2000 iterations, which relates to approximately 5-10 
oscillations of the wake depending on the condition. The solutions are not periodic and are not oscillating 
about a mean value. Although the results will not be shown here, a more detailed study of this model and 
its convergence behavior was done. It was found that the model does not produce the expected periodic 
behavior and the model has difficulty converging to a mean value. This behavior is most likely caused by 
the lower value of a 0 used in the lagRST-2 model ( a 0 = 0.2). For the lagRST implementation, this choice of 
coefficient causes convergence issues and is not recommended for future use. 

II. G. 3. Comparison with Experiment 

Table 14 shows the Cl and Cd results for all turbulence models run on the medium grid and the percent 
difference from the wind tunnel results. The SA and SST models by far have the highest percent error at the 
lower a’s. Figure 28 compares Mach contours for the the Spallart Allmaras and SST results at a = 142.4° 
to the Lag model, which has the lowest percent error. The SA model has a noticeably thicker wake. This 
larger wake increases the disturbance to the flowfield, and thus increases the Cd (indicated by a very large 
negative percent error). The Spallart Allmaras and SST models also both have stronger shocks than the 
Lag model. This also accounts for higher drag by reducing P t . For the SA and SST models, the precent 
error reduces as the a increases. However, for the SA model, even at a = 170.4° the predictions are still off 
by approximately 20%. The baseline SST model provides better results than the Spalart model at all a’s, 
but still has a percent error between 12 - 16 % for the two low a cases. The prediction for a = 170.4° is 
excellent, with percent error around 5%. 

When all three a’s are considered, the Lag model provides the best results. The Cd results are nearly 
identical at the lower a’s, and and the Cl results have some of the lowest percent errors at the same 
conditions. At a = 170.4°, the percent errors raise significantly over the lower a numbers, but still provide 
a prediction less than 10%. 

For the lagRST-1 and lagRST-3 models, the percent errors are of the same order for all a’s. This is a 
definite improvement over the SST model for the lower a cases, but for the a = 170.4° case the errors are 
larger than the SST model. It is also worth noting that the simulations predict a Cd and Cl above the 
experimental results for the lower a cases, but predict a value below the experimental results for the highest 

a. 


28 of 37 


American Institute of Aeronautics and Astronautics 




(a) Lag 


(b) lagRST-1 
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(c) Mach Contour Legend 


Figure 26. Mach contours in the pitch plane for the lag and lagR.ST-1 models showing the wake at a = 142.4°. 
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Figure 27. Time history of Cd (plotted vs iteration) for the lagRST-2 model at all o’s. 
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(a) SA 


(b) SST 



(c) Lag 


(d) Mach Contour Legend 


Figure 28. Mach contours in the pitch plane for the SA and SST models showing the wake at a. = 142.4°. 
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The lagRSTSST has the best prediction of any model at a = 170.4°, and the percent errors at the two 
lower a’s are as good or better than any of the other models, excluding the Lag model. Even though Cd 
and Cl results are very close for each a, the predictions have a non-oscillating wake. The lagRSTSST model 
is a huge improvement over the baseline SST, especially at the lower a’s, because the percent error for Cd 
reduces by almost 10% at a = 170.4°. Although the SST model has a small percent error to begin with, 
the lagRSTSST reduces the % error by around a factor or two over the baseline SST model. Figure 29 
shows four different predictions with SST based models. As previously discussed, the baseline SST provided 
an unsteady wake, whereas the lagRSTSST does not (this is compared in figure 29(a) with figure 29(c)). 
At first glance, this behavior is inconsistent with the assumed behavior of the lag-based models. However, 
there is a valid explanation. The baseline SST model used in these results is the implemented version in 
OVERFLOW 2.2c, however the SST model used to define the equilibrium conditions for the lagRSTSST 
model is taken from OVERFLOW 2.0aa. Figure 29(b) shows an SST result with the OVERFLOW 2.0aa 
implementation, but with the new version of the code. This explains the previous inconsistency. The newer 
implementation of the SST model produces an unsteady answer. The differences are in the form of limiters 
in the 2.2c version of the code that affect the value of cu in the source calculations, production of k, and the 
cross diffusion. 

Because the lagRSTSST results were steady, a second set of cases using a lower value of a 0 value (a 0 = 0.2) 
were run. With the lower value of a 0 , it was postulated that the results could become unsteady and thus 
more representative of the wind tunnel test. The results for this test case are labeled “lagRSTSST low a 0 ” 
in table 14. The Cl and Cd results are nearly identical to the lagRSTSST model. However the force and 
moment data exhibit a periodic oscillation. Figure 29 shows the Mach contours comparing the lagRSTSST 
model with the nominal and lower value of a 0 . Although the Cd and Cl results are very similar, the wake 
contours are not. Figure 29(d) shows that by reducing the value of a 0 , the wake near the point where 
re-circulation would close becomes unsteady and sheds. This is further downstream than what the baseline 
SST model shows and is slightly upstream of what the Lag model predicts. 

Table 15 shows the computed Strouhal numbers for the cases that had oscillating wakes. The number in 
parenthesis is the percent of the amplitude of the oscillation referenced to the Cd- The Strouhal number for 
the numerical simulations was extracted by analyzing the Cd and Cl time varying results. The difference 
in physical time between peak to peak oscillations in Cd was used to determine the frequency of oscillation. 
The results was then verified with Cl , which produced nearly identical results. For the wind tunnel data, 
time varying Cd or Cl results were not available, so the Strouhal number was extracted from unsteady 
pressure data taken from pressure transducers on the backshell. None of the predictions had oscillations 
at 170.4°. The lagRST-1, lagRST-3, Lag, and SST models all had a repeatable oscillation pattern in the 
wake for a’s of 142.4° and 154.4°. The lagRSTSST model did not have an oscillating wake when using an 
a 0 value of 0.35, but when it was reduced to 0.2, the solution oscillated around a mean value. However, 
the magnitude of the oscillation was approximately 0.01% of the total Cd- There was no oscillation at 
a = 154.4° for this model. In general, the results provided by the models are comparable between models. 
At the lowest a, the lagRST-1, lagRST-3, and SST models all have about the same prediction. This is 
not true for the lagRSTSST (using low a 0 ), which predicts a lower value that is more comparable to the 
experimental results. However at the two higher a’s, there is no oscillation. 

III. Summary 

The goal of the current research was to advance current turbulence modeling capabilities in the pre- 
diction of shock wave turbulent boundary layer interactions and flows with massive separation for complex 
configurations, relevant to NASA. The current methodology involved taking a baseline turbulence model 
(k — co or SST) and using it to define an equilibrium turbulent condition. The Reynolds stresses were then 
relaxed towards this equilibrium condition with a simple “lag” equation, thereby mimicking the expected 
behavior through which the streamline history influences the turbulence strain rate. This new turbulence 
model addresses deficiencies observed in flows where varying adverse pressure gradients are seen. The novel 
turbulence model extends work by Olsen and Coakley' and Olsen et al. 19,29 to include Reynolds stress 
transport. By solving for the Reynolds Stresses and not the turbulent eddy viscosity, the models are allowed 
to relax to their non-equilibrium values with more degrees of freedom. This work represents one of the 
first set of computations with a robust Reynolds stresses model in a production CFD code applied to real 
geometries of interest in a URANS method. 


31 of 37 


American Institute of Aeronautics and Astronautics 










•=' BPBB 


BKlfl 

■r 

B 






Bl 







































IB 

P 

■ 




BIB 

BUB 

WA 

Iki 

iJfliiB 

TM 

fl 

B 




B 

Bl 

wmm 

d BB 

aa 

Bl 

aa 

BB 

■ k'^Baa 
BBMBB 

BiBiaaiBir.aiaiaa 
BB BIBrjBBPBBBBBP* 

■ a i B a a a ' jiaa 

BBIllBriBBBB 

■ a 
BB 

aa 

BB 

aaaaaaaa 

IBBBBBBB 

aa 

BB 

BBI 

BBI 


■ 


■ bbbbbbbbbi'bbbb ii a j ■■■■■■■■■■■■■■■■■■■■■■■ 


■ ■■ ■■■■ aar ■ ar.aa 


BBt'a Bl IB ■ a ■ a ■ 


¥ 


W& 


■ 

■ 

K 




ana imnii».«inini»i.ii «u hi.— t i n'l i t 

-■■■■kiaiaaiaraii la ■iBiBaia»-aiii a ■i'aaiB^aiii ■ ■iiBftHi' , i'k.Ba ivi : aa.% 

■ ■■■■■K?«ii>jiBiviaiiaiakiBiaiiriMaiaFiiaaiBii aa iia^k-rr ■ mu — 

■MMBLVBBBr ia BIBB rBBBBI BBB1B1 BKBCZ»<IOr.»' 

■ BaaiaBaaaaaaaiiaia aaiiaa-aaa.'BPa jaa aaia! aaaanaaai. a a.* at < r .fJ'<r 

■ ■■BBBBBBIBfllfll IBB BIB IB^IBIBBIBI BB BBI BBBBI BIBB IBBk'H^drj'ir/ 


IBBBBBBRBBIB.lKRBB'flllB'IBLlBBEP&ffli 
IBBBBBR 'BBI WUBK^BIIBlIBKBBBBRIll 

kaiiipiiiiniiuiip;j|i«iiB» r m w i a 

(a) SST OVERFLOW V2.2 



(b) SST OVERFLOW V2.0 



(c) lagRSTSST 


(d) lagRSTSST low a B 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 


(e) Mach Contour Legend 


Figure 29. Mach contours in the pitch plane for the SST and lagRSTSST model with the nominal and lower values of 
a 0 showing the wake at a — 142.4°. 
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Table 14. Cl and Cd results and percent errors for all turbulence models on the medium grid. 


Model 

Alpha 

C D 

C L 

Cd% error 

Cl% error 

SA 

142.4° 

1.04 

0.71 

-25.8 

-31.5 

SA 

154.4° 

1.23 

0.58 

-23.8 

39.2 

SA 

170.4° 

1.39 

0.23 

-19.4 

-18.2 

SST V2.2 

142.4° 

0.95 

0.63 

-15.2 

-16.2 

SST V2.2 

154.4° 

1.12 

0.49 

-12.6 

-16.1 

SST V2.2 

170.4° 

1.11 

0.18 

4.1 

5.5 

Lag 

142.4° 

0.83 

0.58 

-0.9 

-7.1 

Lag 

154.4° 

0.98 

0.44 

0.7 

-4.4 

Lag 

170.4° 

1.09 

0.17 

6.3 

10.0 

lagRST- 1 

142.4° 

0.87 

0.59 

-5.7 

-9.0 

lagRST- 1 

154.4° 

1.06 

0.46 

-6.7 

-10.2 

lagRST-1 

170.4° 

1.08 

0.17 

6.7 

12.1 

lagRST-3 

142.4° 

0.87 

0.59 

-5.2 

-7.8 

lagRST-3 

154.4° 

1.06 

0.46 

-6.9 

-8.7 

lagRST-3 

170.4° 

1.09 

0.17 

6.2 

12.2 

lagRSTSST 

142.4° 

0.85 

0.59 

-3.1 

-8.6 

lagRSTSST 

154.4° 

1.03 

0.46 

-4.4 

-10.1 

lagRSTSST 

170.4° 

1.19 

0.19 

-2.3 

0.9 

lagRSTSST low a Q 

142.4° 

0.83 

0.57 

-4.5 

-7.4 

lagRSTSST low a Q 

154.4° 

1.01 

0.44 

-2.3 

-6.0 

lagRSTSST low a 0 

170.4° 

1.16 

0.19 

-0.3 

4.1 

SST V2.0 

142.4° 

0.88 

0.62 

-7.1 

-13.8 

SST V2.0 

154.4° 

1.07 

0.49 

-8.3 

-17.3 

SST V2.0 

170.4° 

1.22 

0.20 

-5.6 

-4.8 


The initial study was on four computationally “simple” test cases, an incompressible flat plate, a turbulent 
mixing layer, an incompressible adverse pressure gradient turbulent boundary layer with separation, and a 
transonic bump. The results of this initial study showed mixed results, with no singular implementation of 
the model providing the best results. All of the implementations worked well for the incompressible turbulent 
boundary layer, but none of the models were able to predict the proper shape of the velocity profile near the 
edge. This is not surprising, as typical closure techniques concentrate on matching the wall properties and 
the law of the wall first. For the mixing layer, the k — ui based models incorrectly predicted the spreading 
rate by nearly a factor of two below the experimentally determined value. In contrast, the SST based models 
predicted a spreading rate that matched the experimental data very well. This result was unexpected and 
is still being assessed. In contrast, the separation predictions were far more accurate using the k — uj based 
models, with the standard lagRST model providing excellent results. The SST based models over predicted 
the separation zone for both the incompressible and transonic test cases by over 40%. 

The major set of analysis done in this project centered around the predictions on two test cases. The first 


Table 15. Computed Strouhal number for the lagRST- 1, lagRST-3, SST, and lagRSTSST (with a Q = 0.2) models. The 
numbers in parenthesis after the Strouhal number is a measure of the oscillation. It is the percent the amplitude of 
the oscillation is of the total Cd • 


a 

Data 

lagRST-1 

lagRST-3 

SST 

lagRSTSST 

Lag 

142.4° 

0.22 

0.27 (1.7%) 

0.25 (1.4%) 

0.26 (2.5%) 

0.22 (0.01%) 

0.27 (0.64%) 

154.4° 

0.17 

0.18 (1.8%) 

0.18 (1.9%) 

0.19 (2.1%) 

n/a 

0.24 (0.24%) 

170.4° 

0.15 

n/a 

n/a 

n/a 

n/a 

n/a 
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test case was a hypersonic SWTBLI of an impinging oblique shock on a fully developed turbulent boundary 
layer that caused separation. For this test case, the baseline SST model provided separation zone length 
predictions that were less than 1% from the experimentally observed values. The lagRSTSST and lagRST 
models both predict separation zone lengths on the order of 50% too large. In addition to separation zone 
accuracy, of particular interest too was the peak heat flux prediction at reattachment. The SST model over 
predicted this heat flux by 65%. By lagging the Reynolds stresses, the peak heat flux prediction reduced 
by 25% for the lagRSTSST. So although the lagging technique had an adverse affect on the prediction of 
the separation zone, it did have a positive effect on one of the negative characteristics of the SST model for 
SWTBLI. 

The final analysis was done on the Orion (now MPCV) capsule. The grid system and techniques used here 
are the ones used by NASA’s aerodynamics teams building the aerodynamic databases for these vehicles, and 
the baseline models are the ones used by NASA. The first contribution of this analysis is that the Reynolds 
stress modeling implemented in the OVERFLOW code is robust and can be run in an unsteady, time accurate 
fashion and provide grid converged results. This is a significant contribution since the predictions are solving 
for Reynolds stresses. The results for this case showed excellent results for the lag technique. The standard 
turbulence models (Spalart Allmaras and SST) almost always provided higher percent errors than the lag 
based models. The Lag, lagRST, and lagRSTSST models all were able to generally predict the Cn and 
Cl to less than 10%. Several of the lag based models predicted oscillating wakes with Strouhal numbers 
comparable to the experimentally obtained values at lower a’s. Unfortunately, none of the models were able 
to predict an unsteady wake at the highest a. 

IV. Recommendations 

It is clear from the analysis done that a simple lag based Reynolds Stress model does not solve all the 
problems identified for flowfields with massive separation or SWTBLI interactions. Based the suite of cases, 
two models are recommended for further study. The lagRST-1 model should be treated as the baseline 
lagRST model, with the same coefficients used in the standard eddy viscosity Lag model, a Q = 0.35 and 
CTfc = 0.8. This model did very well for incompressible and transonic flow predictions. It should be noted that 
if the flow of interest has a mixing layer component, the predictions need to be validated with experiments 
because of its poor predictions of the growth rate of the Bell Mehta mixing layer. Although this mixing layer 
issue needs to be explained, it did not seem to affect the capsule predictions. The lagRSTSST model with 
the standard coefficients ( a 0 = 0.35) also performed well for the complicated test cases. The issue with this 
model is its separation prediction for incompressible and transonic flat plates, where the results were very 
poor. The separation predictions for this model were actually better for the SWTBLI than for the simpler 
cases, and the lag equations reduced the peak heat flux at the reflected shock. These are promising results, 
as an over prediction of the magnitude SST had would most likely result in over design of the TPS system. 

The SWTBLI results should be verified with other wind tunnel datasets to ensure consistent behavior. 
Brown 34 has identified two other high quality experimental datasets that could be used. One dataset is 
the Mach 8.9 experiment of Murray 46 ™ 48 which was conducted at the Imperial College Nitrogen gun tunnel. 
This test series consists of a hollow axisymmetric cylinder with an axisymmetric cowl use as a shock ring 
generator. Data available for this case are wall pressure and heat flux. The second experiment is the Mach 
8.18 experiment of Kussoy and Horstman. 49 It was conducted on a 2-D flat plate in the no longer operational 
NASA Ames Research Center Hypersonic facility. The test configuration was the nominal 2-D flat plate 
with a shock generator inclined, similar to the Schulein configuration assessed in this study. Pressure, skin 
friction, and heat flux were obtained in this study as well as flow field quantities (no turbulence quantities) 
were obtained. In addition, these configurations should be assessed for 3-D affects. 

In addition to simulating more configurations, there are also several other forms of the turbulence model 
that should be modeled. In particular, the form of the Reynolds stress tensor used in the lagRST model 
could provide numerical and / or physical modeling improvements over the current implementation. Two of 
these forms are detailed below. The anisotropy tensor has been used by several other researchers 22,23 ’ 50 ’ 51 
and may provide improvements to the predictions. 

Instead of using r,:j , the mean normal stress component can be subtracted from the Reynolds stress to 
leave the deviatoric portion of the tensor. The Boussinesq equation ( 6) can be rewritten to show 
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where the superscript D denotes the cleviatoric part of the tensor. The lag equation would then become the 
following. 
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In this formulation, the turbulent kinetic energy solved for in equation 3 would be the actual lagged variable 
used in all flowfield computations. 

The Reynolds stress anisotropy tensor is another option for the variable to be lagged. It is defined here. 
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By substituting equations 2 and the relation Tkk = ~P u 'k u 'k = ~2 pk into the equation for bij, you arrive at 
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The numerator of the previous equation is the deviatoric Reynolds stress tensor (left hand side of equation 13), 
thus substitution leaves 
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Equation 14 can then be used to define the bij e , the equilibrium anisotropic Reynolds stress tensor, 
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Another possibility is to remove the Boussinesq approximation and drive the Reynolds Stress tensor to 
a higher closure model. This could be done by solving for an algebraic Reynolds Stress tensor. 

Coakley 52 proposed what he refers to as the KEA (the letters KE indicates that this version of his model 
tries to follow k — e model behavior and the A represents the choice of switching function used) model, 
which was intended to address deficiencies in the SST model. Coakley’s primary purpose was to remove the 
dependence on wall distance. This would improve the timings of the lag based models based on the SST 
formulations, while still keeping the accurate mixing layer predictions that eluded the k — to based models. 
Coakley shows results that are extremely close to the SST predictions with the KEA model. 
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